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Abstract: The main purpose of this paper is to provide a method of multi-criteria decision making that 
combines interval neutrosophic sets and TOPSIS involving the relative likelihood-based comparison 
relations of the performances of alternatives that are aggregated into interval numbers. A fuzzy cross-entropy 
approach is proposed to state the discrimination measure between alternatives and the absolute ideal 
solutions after a transformation operator has been developed to convert interval neutrosophic numbers to 
simplified neutrosophic numbers. Finally, an illustrative example is given and a comparison analysis is 
conducted between the proposed approach and other existing methods to verify the feasibility and 
effectiveness of the developed approach. 
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1. Introduction 

In practice, multi-criteria decision making (MCDM) methods are widely used to rank alternatives or 
select the optimal one with respect to several concerned criteria. However, in some cases, it is difficult for 
decision-makers to explicitly express preference in solving MCDM problems with uncertain or incomplete 
information. Under these circumstances, fuzzy sets (FSs), proposed by Zadeh [1], where each element has a 
membership degree represented by a real number in the standard interval [0,1] , are regarded as a significant 
tool for solving MCDM problems [2,3], fuzzy logic and approximate reasoning [4]. Sometimes, FSs cannot 
handle the cases where the membership degree is uncertain and hard to be defined by a crisp value, therefore, 
the concept of interval-valued fuzzy sets (IVFSs) was proposed [5] to capture the uncertainty of membership 
degree. Generally, if the membership degree is defined, then the non-membership degree can be calculated 
by default. In order to deal with the uncertainty of non-membership degree, Atanassov [6] introduced 
intuitionistic fuzzy sets (IFSs) which is an extension of Zadeh’s FSs, and the corresponding intuitionistic 
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fuzzy logic [7] was proposed. IFSs consider both the membership degree and the non-membership degree 
simultaneously. So IFSs and intuitionistic fuzzy logic are more flexible in handling information containing 
uncertainty and incompleteness than traditional FSs. Currently, IFSs have been widely applied in solving 
MCDM problems [8- 16]. Moreover, in reality, the degree of membership and non-membership in IFSs 
may be expressed as interval numbers instead of specific numbers. Hence, interval-valued intuitionistic fuzzy 
sets (IVIFSs) [17] were proposed, which is an extension of FSs and IVFSs. In recent years, MCDM 
problems with evaluation information derived from IVFSs have attracted much attention of researchers 
[18- 24], in which averaging operators, aggregation operators, prospect score function and possibility 
degree method are involved. Furthermore, TOPSIS, proposed by Hwang and Yoon [25], was also used for 
MCDM problems under IVIFSs environment in [26,27]. 

Although FSs and IFSs theory have been developed and generalized, they cannot deal with all sorts of 
uncertainties in real problems. Some types of uncertainties, such as the indeterminate and inconsistent 
information, cannot be handled. For example, a paper is sent to a reviewer, he or she says it is 70% 
acceptable, 60% unacceptable and the statement is 20% uncertain. This issue cannot be handled effectively 
with FSs and IFSs. Therefore, some new theories are required. 

Smarandache [28,29] coined the concept of neutrosophic logic and neutrosophic sets (NSs). Rivieccio 
[30] pointed out that a NS is a set where each element of the universe has a truth-membership, 
indeterminacy-membership and falsity-membership, respectively, and it lies in the non-standard unit interval 

]CT, l + [ . The uncertainty presented here, i.e. indeterminacy-membership, is quantified and independent of 

truth-membership and falsity-membership, which is different from the argument that the incorporated 
uncertainty is dependent on the degree of membership and non-membership of IFSs [3 1 ]. Thus, it is practical 
and flexible to address problems containing uncertain, incomplete and inconsistent information with NSs. 
For the aforementioned example, it can be presented as .*(0.7, 0.2, 0.6) by means of NSs. In recent years, 

the NS theory has found practical applications in various fields, such as relational database systems, semantic 
web services [32], mineral prospectivity prediction [33], image processing [34-36], granular computing [37], 
medical diagnosis [38] and information fusion [39]. It is clear that all of those proposals, promising as they 
are, still need to be refined from a formal point of view. Without specific description, it is hard to apply NSs 
in real scientific and engineering situations. Hence, the notion of single-valued neutrosophic sets (SVNSs) 
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was introduced by Wang et al. [40], which is an instance of NSs. Majumdar et al. [31] introduced similarity 
and entropy measures of SVNSs. Furthermore, the correlation and correlation coefficient of SVNSs [41] as 
well as the cross-entropy of SVNSs [42] were presented. In addition, Ye [43] also introduced the concept of 
simplified neutrosophic sets (SNSs), which can be described by three real numbers in the real unit interval 
[0,1]. Peng et al. [44] defined the novel operations and aggregation operators of SNSs and applied them in 
solving multi-criteria group decision making problems. 

Actually, the degree of truth, falsity and indeterminacy about a certain statement are denoted by several 
possible interval numbers instead of three real numbers in SNSs in some real situations. For example, an 
expert is asked to give the opinion about a certain statement, he or she may provide his or her own 
evaluations based on surveys as well as his or her knowledge and experience, and the evaluations are 
gathered in the form of interval numbers. The possibility that the statement is true is between 0.5 and 0.7 and 
that the statement is false is between 0.2 and 0.4 and the degree that he or she is not sure is between 0. 1 and 
0.3. That is beyond the scope of SNSs. So the concept of interval neutrosophic sets (INSs) was proposed by 
Wang et al. [32], which is a particular instance of NSs. For the aforementioned example, it can be presented 

as x([0. 5,0. 7], [0. 1,0.3], [0.2, 0.4]) by means of INSs. Subsequently, the studies of INSs had been 

conducted in various aspects, which concentrated mainly on defining and amending operational laws [45,46], 
correlation coefficients [47], similarity measures [48] and distances [49] to aggregate opinions of experts or 
decision-makers in MCDM problems, where the opinions of experts are expressed by INSs. Nevertheless, it 
is complex to directly process the assessment information with INSs. Furthermore, most of current studies 
coping with MCDM problems with INSs by means of aggregating interval numbers into real numbers, which 
may bring some drawbacks in operations and lead to great information loss [27,50]. 

As an important topic in the theory of FSs, entropy measure of FSs mentioned by Zadeh [1], is an useful 
tool to measure uncertain information. Burillo and Bustince [51] introduced the notion that entropy of IFSs 
and IVIFSs can be used to assess the degree of intuition of an IFS or IVIFS. In the course of determining 
criterion weights under intuitionistic fuzzy environment, Chen and Li [52] utilized the intuitionistic fuzzy 
entropy measure to estimate the objective criteria weights. Despite the existing research effort, dealing with 
MCDM problems with completely u nkn own criteria weights in the framework of IVIFS remains an open 
problem [53,54]. The starting point for the cross-entropy method is information theory as developed by 
Shannon [55]. Cross-entropy is applied to measure the discrimination information, according to Shannon’s 
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inequality [56]. Then Shang and Jiang [57] proposed a fuzzy cross-entropy and a symmetric discrimination 
information measure between two FSs. Furthermore, the concepts of intuitionistic fuzzy cross-entropy and 
discrimination for IFSs were introduced by Vlachos and Sergiadis [58]. A fuzzy cross-entropy on 
interval-valued intuitionistic fuzzy numbers was proposed by analogy with the intuitionistic fuzzy 
cross-entropy [59]. A series of mathematical programming models based on cross-entropy were constructed 
to determine the weights of criteria under IVIFSs environment [27]. Other applications with the fuzzy 
cross-entropy include portfolio selection [60], the divergence of uncertain variables [61], MCDM problems 
[62-64] and so on. 

Flowever, compared with fuzzy cross-entropy that is widely applied in FSs and IFSs, there exists few 
studies on the application of entropy theory under interval neutrosophic environment and attention paid to 
avoid information loss is not enough in the process of information aggregation. What’s more, TOPSIS which 
is one of popular decision making methods has been applied effectively and availably in practice to address 
problems in various fields. For this purpose, a MCDM method on the basis of cross-entropy is proposed, in 
which a transformation operator is defined to convert each interval neutrosophic number (INN) to a 
simplified neutrosophic number (SNN). It can effectively avoid artificially setting parameter values of the 
transformation operator in IVIFSs and IFSs [54,65] and to some degree, lower the computational complexity 
in directly processing evaluation information with INSs at the same time. In addition, a TOPSIS method 
associated with possibility degree method is proposed to rank the performances of alternatives in the form of 
interval numbers, which can decrease the loss of assessment information and guarantee the reasonability of 
the ultimate ranking results. 

The rest of the paper is organized as follows. In Section 2, interval numbers as well as the concepts and 
operations of NSs, SNSs and INSs are briefly reviewed. In Section 3, an operator is put forward, which can 
transform each INN into a SNN, together with a cross-entropy for SNNs and some useful properties. 
Subsequently, TOPSIS for MCDM problems with INSs is developed in Section 4. In Section 5, an 
illustrative example is given and a comparison analysis is conducted between the proposed approach and 
other existing methods to verify the feasibility and effectiveness of the developed approach. Some summary 
remarks are given in Section 6. 
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2. Preliminaries 



In this section, some basic concepts and definitions related to INSs, including interval numbers, 
definitions and operational laws of NSs, SNSs and INSs are introduced, which will be utilized in the latter 
analysis. 

2.1 Interval numbers 

Interval numbers and their operations are of utmost significance to explore the operations for INSs. In 
the following, some definitions and operational laws of interval numbers are given. 

Definition 1 [66,67]. Let a = [a L ,a u ]- {x\a L <x<a u } , then a is called an interval number. Especially, 
a - [a L , a u ] will be degenerated to a real number if ci L = ci u . 

Consider any two nonnegative interval numbers a = [a L ,a u \ and b=\b L ,b u ], where 0 <a L <x<a u , 
0 < b L < x < b u . Then their operations are defined as follows [66]: 

(1 )a=boa L -b L ,a u =b u ; 

(2) a + b — [a L + b L , a u + b u ]■ 

(3 )a-b=[a L -b u ,a u - b L ]; 

(4) Aa — \Aa L ,Aa u ],A > 0. 

Definition 2 [68,69]. For any two interval numbers a-[a L ,a u ] and b -\b L ,b u ], the possibility of a>b 

is formulated by p{a >b) = max] 1 -maxi — Q _ ,0 1,0 > , where L(a) = a u -a L and L(b)-b u -b L . 

{ { L(a) + L(b ) J J 

The possibility degree of a > b has the following properties [68]: 

(1) 0 < p(a>b) < 1; 

(2) p(a > b) - p(b > a) = 0.5 if p(a>b)=p(b>a ); 

(3) p(a >b) + p(b > a) = 1; 

(4) p(a>b) = 0 if a u <b L , p(a>b) = 1 if a L >b u ; 
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(5) For any interval numbers d, b and c, p(d>c)> 0.5 if p(a>b)> 0.5 and p(b>c)> 0.5; 
p(d >c) = 0.5 if and only if p(a >b)= p(b >c) = 0.5 . 



2.2 NSs and SNSs 

Definition 3 [28], Let X be a space of points (objects) with a generic element in X , denoted by x . A NS 
A in X is characterized by a truth-membership function t A (x) , an indeterminacy-membership function 

i A (x) and a falsity-membership function f A (x) . t A (x) , i A (x) and /, (x) are real standard or 

nonstandard subsets of ]0“, l + [ , that is, t A (x) : X — »]0“, l + [ , i A (x) : X — >]0“, l + [ , and 

f A (x) : X — »]0“, l + [ . There is no restriction on the sum of t A (x) , i A (x) and / A (x) , so 

0 < sup t A (x) + sup i A (x) + sup /, (x) < 3 + . 

Definition 4 [28], A NS A is contained in the other NS B, denoted by AczB if and only if 
inf t A (x)< inf t B (x) , sup t A (x) < sup t B (x) , inf i A (x) > inf i B (x) , sup i A (x) > sup i B (x) , inf /, (x) > inf f B (x) 

and sup f A (x) > sup f B (x) , for any x e X . 

Since it is hard to apply NSs to practical problems, Ye [43] reduced NSs of nonstandard interval 
numbers into a kind of SNSs of standard interval numbers. 

Definition 5 [30,43], Let X be a space of points (objects) with a generic element in X , denoted by x. A 
NS A in X is characterized by t A (x ) , i A (x) and f A (x ) , which are single subintervals/subsets in the 

real standard [0,1], that is, t A (x) : X — > [0,1] , i A (x ) : X —>[0,1], and f A (x) : X — » [0,1] . And the sum of 

t A (x) , i A (x) and f A (x) satisfies the condition 0 < t A (x) + i A (x) + f A (x) < 3 . Then a simplification of A 

is denoted by A = { ( x, f ^ ( x), i A ( x), f A (x) ) | x el}, which is called a SNS. It is a subclass of NSs. If ||x||=l , a 

SNS will be degenerated to a SNN, denoted by A = (t A ,i A ,f A ) . 

Definition 6 [30,43]. A SNS A is contained in the other SNS B, denoted by AcB if and only if 

t A (x) < t B (x) , i A (x) > i B (x) and f A (x) > f B (x) , for any xeX . Especially, A = B if Acfi and BcA. 
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The complement set of A, denoted by A c is defined as A c ={(x,f A (x),i A (x),t A (x))\x <= A} . If A is a 
SNN, then A c ={(f A ,i A ,t A )) . 



2.3 INSs 

In the actual applications, sometimes, it is not easy to express the truth-membership, 
indeterminacy-membership and falsity-membership by crisp values, and they may be easily described by 
interval numbers. Wang et al. [32] further defined INSs. 

Definition 7 [30,32], Let X be a space of points (objects) with generic elements in X , denoted by x. An 
INS A in I is characterized by a truth- membership function t- A (x) , an indeterminacy-membership 

function (x) and a falsity-membership function /Ax). For each point x in X , we have that 



t A {x) = [t L Ax),t u ,{x)\ 






/iW = [//(4/f(x)]c[0,i] 



and 



0 <t A (x) + / - (x) + /- (x) < 3 . We only consider the subunitary interval of [0,1] . An INS is the subclass of a 
NS. Therefore, all INSs are clearly NSs. Especially, an INS will be reduced to a SNS if t A (x)-t A (x), 
i\{x) = i' A (x) and f A (x) = f A (x) . In addition, if |x|=l , an INS will be degenerated to an INN, denoted by 

Definition 8 [30,32]. An INS A is contained in the other INS B , denoted by AcB if and only if 
t L x {x)<t L ,{x), A y (x)<^(x), ; '(x)>^(x), i u A (x) > (x) , // (x) > // (x) and //(x)>//(x), for any 

xeX. Especially, A- B if AcB and BcA. The complement set of A , denoted by A c is defined as 
A c ={(x,[//(x),//(x)],[/|(x),/J(x)],[^(x),^(x)])|xeZ} .If A is an INN, then 

a c ={([//, 

Definition 9 [46]. Let A = (\t L k ,t u k ],[i L A ,i u x ],[//, f u A ]) and B = ([^MM ]> [//.//]) two INNs, 
where A > 0 . The operations for INNs are defined as below. 
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(2) A' =([((‘)',(r? )'],[! -(!-/" ) J ],[1-(1 



(3)A + fl = ([<i+rJ-<P<J,<J+<i[-. 



■<a AuTM <wi fi.fj! -f^y. 

(4) a ■ B = (lit ■ 4 , r» ■ ],pf + i‘ - ;f ■ /{ , i? + - i" . /■' 

Example 1. Assume two INNs A = ([0.7, 0.8], [0,0.1], [0.1, 0.2]) and B - ([0.4, 0.5], [0.2, 0.3], [0.3, 0.4]) , and 
X = 2 . Then the following results can be calculated. 

(1) 2- A = ([0.91, 0.96], [0,0.01], [0.01, 0.04]); 

(2) A 2 = ([0.49, 0.64], [0,0.19], [0.19, 0.36]); 

(3) A + B = ([0.82, 0.90], [0,0.05], [0.03, 0.08]); 



(4) A - B = ([0.28, 0.40], [0.20, 0.37], [0.37, 0.52]). 



Bustince and Burillo [70] put forward an operator H p , which can transform each IVFS into an IFS. 
As an improved extension of H p , an operator H p r is defined to convert each INN to a SNN. 

Definition 10. Let p,q,r e [0, l] be three fixed numbers, any INN can be transformed into a SNN through 
the operator H p q r . 



H 



M) = (i l i + p-KA + ^ w >,’fi + r - w ,,)’ 



where W t =i A -t- , W t - i. - i, and W f =/- -/- . Obviously, H (A) is a SNN and can be 



determined with respect to p , q and r . That is, H p r (A) is well defined in all value ranges of p , q 
and r . 

Example 2. Use the data of Example 1. Then the following transformation results can be obtained utilizing 
the operator H p r . 



(1) H p q r {A) = (0.7 + p ■ 0.1,0 + q ■ 0.1, 0.1 + r • 0.1); 

(2) H p q r (B) = (0.4 + p ■ 0. 1,0.2 + q ■ 0. 1,0.3 + r ■ 0. l) 
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It is shown that H (A) and H pqr {B) will be two specific SNNs if the values of p , q and r 
are given. And the method to determine the parameter values will be discussed in detail in Section 4. 

3. Cross-entropy for SNNs 

In an analogous manner to the proposals of Ye [42] and Vlachos and Sergiadis [58], the following 
definition of fuzzy cross-entropy for SNNs is proposed. 

Definition 11. Let A and B be two SNNs. Then the cross-entropy between A and B can be defined as: 

I NS (A,B) = t A ln-^- + i A In— — 7 - + f A In . (1) 

t A +t B l A +l B + /fl 

Eq. (1) can indicate the degree of discrimination of A from B . It is obvious that I NS (A,B) is not 
symmetric with respect to its arguments. Therefore, a modified symmetric discrimination information 
measure based on I NS (A, B) can be defined as: 

D ns (A,B) = I ns (A,B) + I ns (B, A) . (2) 

The larger D NS (A,B ) is, the larger the difference between A and B will be, and vice versa. 
Property 1. Let A and B be two SNNs. Define the degree of discrimination of A from B as 
D ns (A,B) . Then the following properties hold. 

(1 )D ns (A,B) = D ns (B,A), 

(2) D ns (A,B) = D ns (A c ,B c ), where A c and B‘ are the complement sets of A and B , respectively, 
defined in Definition 6 . 

(3) D ns (A, B) > 0 and D NS (A, B) = 0 if and only if A=B . 

Proof. Obviously, it can be easily verified that (1) and (2) hold. In the following, the proof of (3) is shown as 
below. 

Now consider the function f(x) = x In x, where v e (0,ll . Then f'(x) - 1 + lnx and f"(x) = — >0, 

x 

where xe(0,l] . Accordingly, f(x) = x lnx is a convex function. Therefore, for any two points 
x,,x 2 e ( 0 ,l] , the inequality > /( ' Y ' + Xl - ) holds. 
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X ~\~ X 

Utilize f(x) - xlnx in the above inequality and x, lnxj + x 2 lnx 2 - (x, + x 2 ) In ’ 1 2 > 0 can be got, 



where the equality holds only if x, = x 2 . Similarly, the following equation can be obtained. 



t > +t a 



^ns (A,B) — I NS (A, B) + I NS {B, A) — (t A lnf A + t B lnf B ) (f A +t B )ln + 



(i A In i A + i B ln/ B ) - (i A + i B ) ln^liiL + 



(/>A + / B ln/ fi )-(/ A + / B )li 



n 



/a +/b 



Because T> 0, 7>0 and F> 0, D NS (A,B)> 0 holds. D NS (A,B) = 0 holds only if t A -t B , 
i A = i B and f A = f B , namely A = B . 

Example 3. Assume two SNNs A = (0.8, 0.1, 0.2) and B = (0.5, 0.3, 0.4) . Then the following result can be 
obtained by applying Eqs. (1) and (2). 

D ns (A,B) = I NS (A,B) + I ns (B, A)=0. 1212 . 



4. MCDM approach based on cross-entropy and TOPSIS 



This section presents an approach for MCDM problems with INSs by means of the cross-entropy and 
TOPSIS. 

For a MCDM problem, let A = {a l ,a 2 ,---,a m } be a set consisted of m alternatives and let 
C = {cj,c 2 ,-",c n } be a set consisting of n criteria. Assume that the weight of the criterion 



: j (j = 1,2, •••,«), given by the decision-maker, is w., where w ; e [0,1] and =1. And the weight 



;=i 



vector of 



criteria 



can be expressed 



as 



w 



= (w l ,w 2 ,---,w n ) 



Let 



B = 



b;; 



([«].R,?].[C.C])' 



be the decision matrix 



ix, where [r‘,r“], and 



[ft’fX] are interval numbers and represent the degrees of truth-membership, indeterminacy-membership 
and falsity-membership of each alternative a t (z = 1, 2, • ■ • , m) on the criterion cv with respect to the 
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concept “excellent”, respectively. 

In general, there are two types of criterion called maximizing criteria and minimizing criteria. In order 
to uniform criterion types, minimizing criteria need to be transformed into maximizing criteria. Suppose 
the standardized matrix is expressed as R = r. . The original decision matrix B - b,. can be 

L J Jmxn L J Jmxn 

converted to R - r based on the primary transformation principle of Ref. [71 ], where 

L J Jmxn 



hj 



h = ],[iy,iy y_fn-.fi ]), for maximizing criterion c. 

b-j = ])> for minimizing criterion c. 



(3) 



in which /;' is the complement set of b u , defined in Definition 8. 

The absolute positive ideal solution (PIS) and the absolute negative ideal solution (NIS) of INSs 
respectively denoted by a + and a~ , and can be expressed as follows [49]: 

In order to get the cross-entropy or degree of discrimination between cl (i = 1,2, •••,/«) and the ideal 
solutions, each INN is transformed into a SNN based on the operator H p r . Let p tj , q tj , r. e [0, l] be three 
crisp numbers and any INN, denoted by ([thy ],[fy,fy ]) can be transformed into the following 

form: 

H ^.r ( f )=( t « * Pu ■ y . hi + % ■ y , y +r v -w fv ), ( 4 ) 

where U) = , W. = i" - /' and W fii = f" - f u L 



fij J ij J V 



Using Eqs. (1), (2) and (4), the degree of discrimination of (i = l,2,---,m) from a + or a- on 



Cj ( j = 1,2 can be obtained as follows: 



2(th + p„ -W ) 

If -y + P, •" )ln — ^ 1 In- 



y+ P ,-w.+ 1) y+ P , r w, +D 



(/' * r/ -U- )ln2i r U , ) In 2 



(5) 
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2 



D u = (t, + P, • W. ) In 2 + (i + q ■ w. ) In - 






- + ln- 



(tf + V W / f ) ^ 



y '“ J a! i +<l i fW i: +') (< 1 ‘I, • W i, -i) 

2 (tf + v^) 



+ 



(6) 



- + In — , 

+i) (ti +r ir w, +D 



V ‘j ftj 



Using Eqs. (5) and (6), Z> (i can be obtained, which expresses the performance of the 

alternative <3 ■ (i = 1,2,- • -,m) . In other words, the larger D i is, the better the alternative a t will be. 



I) 



ij 



A=y wD =y w, ■ 

^ ' f rr ; /u + /u 



;=i 



(7) 



where D, = — 

7 D- + D,: 



denotes the performance of the alternative a, (/ = 1,2, •••,/«) on c. (j = 1,2, 



n) 



and Wj ( j = 1,2, ••■,«) represents the weight of c . . 

In view of the fact that an INS is characterized by a truth-membership function, a 
indeterminacy-membership function and a falsity-membership function, whose values are intervals rather 
than specific numbers. It is unreasonable to designate a SNN for the given INN by artificially choose only a 
certain p tj , q tj and r.. in Eq. (4) to indicate the evaluation information. Because it may lead to distortion 

or loss of original information [27], Thus, /> may take different values as numerical variations p iJ , q r 

and r tj change. In order to avoid information loss, an interval /> = ^ [)’■ , D'j J is applied to represent the 

performance of the alternative a i (i = \,2,---,m) on c. (j = l,2,---,n) , where and are the 

lower and upper bounds, respectively. 

The values of Z> / ' and // are determined as below. 

As it is shown in Eqs. (5) and (6), both Df and D7. are multivariate continuous functions with 
respect to p tj , q tj and r .. , respectively. They also can reach the maximum and minimum in the domain 

Pii , Qij , r ij e [0, l] • Calculate the partial derivative of D7 and D7. with respect to p tj , q tj and r . , 
respectively. 
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dPa 






Since 2(4 + PiJ ■ W tf ) < t\ + Pij - W ti , + 1 , ^<0. 



ao; sd; 

— 2- = W In 2 > 0 , and — ^ = W, ln2>0. 

8q u j dr * 



5D: 



OD, ; 



Similarly, — ^ = W, ln2>0, — ^ = W, ln- 



2(/,: . ^ u ) 



5/7, 






■(^+^-W +1) 



< 0 . 



and “f -w, ■. y+w 
^ (./;, +0 



< 0 . 



1 ij y fij 

The analysis above indicates that D( can reach its maximum and D7 reaches the minimum if 



p u = 0 and q u = r. = 1 . Likewise, D*. reaches its minimum and D.. reaches the maximum if p u = 1 



and q tj = r = 0 . As a result, it easy to understand that /> = 



D 



d- j+ d; 



can have a minimum when p tj — 0 



and q tj = r.. = 1 and a maximum when p tJ = 1 and q tj - r y = 0 . Then the lower and upper bounds of /> 
can be obtained, respectively. 
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Based on the analysis above, Eq. (7) can be rewritten using operational laws of interval numbers given 
in Definition 1. 



A=2>A = 2>i[ D fD,n 

7=1 j = 1 



2> J D f2>,A 



7=1 



7=1 



( 10 ) 



D j >D j means that a i (i = l,2,---,m) is not inferior to a j (j = \,2,---,m). Then the likelihood 
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matrix P can constructed as below. 



P = (Pij\ 



p 11 


P 12 ' 


■' Pun 


Pll 


Pl2 ■ 


Plm 


Pml 


Pml ■ 


P mm 



(ID 



where p tJ = p(a j > a j ) = p(D j > D j ) = max \ 1 - 



, ( D^-Dl 

1 - max j — = — , 



0 ,0 Land L(D.) - -Df 



(L(D,.) + L(D.) 

According to [69], the ranking vector of the likelihood matrix can be defined as follows: 



co. 



Z m , 
p < + 2 
(m - l) 






, i = 1, 2, ■ ■ ■ , m . 



m 



( 12 ) 



Consequently, the ranking of all alternatives is determined according to the descending order of 
(O' (/ = 1,2, •••,«;) . That is, the larger is, the better the alternative a : will be. 

In a word, the main procedure of the above MCDM approach is listed as below. 

Step 1. Establish the original decision matrix B - /; according to the evaluation information given 

L J J mxn 

in the form of INNs. Then use Eq. (3) to transform B - b r into R - 1 r. I 

L l J J mxn L ” J mxn 

Step 2. Use Eq. (8) and (9) to derive D’: and />' and establish the matrix ( ) = [ D'. , D(' 1 

J J \ J / mxn L J J Jmxn 



Step 3. Use Eq. (10) to calculate the performance of each alternative D j (i = 1, 2, - - - ,m) . 

Step 4. Construct the likelihood matrix P by using Eq. (11) and obtain the ranking vector 
co = (co l ,a> 2 ,---,co m ) based onEq. (12). 



Step 5. Determine the ranking of all alternatives according to the descending order of co t (i = 1, 2, • • • , m) . 



5. Illustrative example 

In this section, an investment appraisal project is used to demonstrate the application of proposed 
decision making approach, as well as the validity and effectiveness of the proposed approach. 

The example on the investment appraisal project of a company adapted from Ref. [48] is employed here. 
ABC company intends to invest a sum of money in the best option. There is a panel with four possible 
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alternatives to invest the money: (1) a l is a car company; (2) a, is a food company; (3) a i is a computer 



company; (4) a 4 is an arms company. This company must make a decision according to the three criteria: 
(1) q is the risk; (2) c 2 is the growth prospects; (3) c 3 is the environmental impact, where c, and c 2 
are maximizing criteria, and c 3 is a minimizing criterion. The weight vector of the criteria is 
w = (0.35,0.25,0.4) . The four possible alternatives are to be assessed under three criteria using INSs, and the 
evaluation information is shown in the following interval neutrosophic decision matrix B — b u 

L J J mxn 

([0.4, 0.5], [0.2, 0.3], [0.3, 0.4]) ([0.4, 0.6], [0. 1,0.3], [0.2, 0.4]) ([0.7, 0.9], [0.2, 0.3], [0.4, 0.5]) 

([0.6, 0.7], [0.1, 0.2], [0.2, 0.3]) ([0.6, 0.7], [0. 1,0.2], [0.2, 0.3]) ([0.3, 0.6], [0.3, 0.5], [0.8, 0.9]) 

([0.3, 0.6], [0.2, 0.3], [0.3, 0.4]) ([0.5, 0.6], [0.2, 0.3], [0.3, 0.4]) ([0.4, 0.5], [0.2, 0.4], [0.7, 0.9]) 

([0.7, 0.8], [0.0, 0.1], [0.1, 0.2]) ([0.6, 0.7], [0.1, 0.2], [0.1, 0.3]) ([0.6, 0.7], [0.3, 0.4], [0.8, 0.9]) 

5.1 Procedures of decision making based on INSs 

In the following, the main procedures of obtaining the optimal ranking of alternatives are presented. 

Step 1. Use Eq. (3) to transform minimizing the criterion c 3 into a maximizing criterion and obtain 

L J _l mxn 

([0.4, 0.5], [0.2, 0.3], [0.3, 0.4]) ([0.4, 0.6], [0. 1,0.3], [0.2, 0.4]) ([0.4, 0.5], [0.2, 0.3], [0.7, 0.9]) 

([0.6, 0.7], [0.1, 0.2], [0.2, 0.3]) ([0.6, 0.7], [0. 1,0.2], [0.2, 0.3]) ([0.8, 0.9], [0.3, 0.5], [0.3, 0.6]) 

([0.3, 0.6], [0.2, 0.3], [0.3, 0.4]) ([0.5, 0.6], [0.2, 0.3], [0.3, 0.4]) ([0.7, 0.9], [0.2, 0.4], [0.4, 0.5]) 

([0.7, 0.8], [0.0, 0.1], [0.1, 0.2]) ([0.6, 0.7], [0.1, 0.2], [0.1, 0.3]) ([0.8, 0.9], [0.3, 0.4], [0.6, 0.7]) 

Step 2. Use Eqs. (8) and (9) to derive D': and />' and establish the matrix ( A, ) 4 3 = [ Ay ] 4 3 • 

([0.4963,0.6597]) ([0.4963,0.8144]) ([0.3317,0.4838]) 

(£>) =\d l D"1 ([° f>954,0.8369]) ([0.6952,0.8369]) ([0.4714,0.7094]) 

l vL 3 L «’ V J 4 x3 ([0.4410,0.6952]) ([0.5433,0.6952]) ([0.5194,0.7146]) 

([0.8369,0.9542]) ([0.6952,0.8902]) ([0.4800,0.5823]) 

Step 3. Use Eq. (10) to calculate the performance of each alternative D i (i = 1,2, 3, 4) . 
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D l = [0.4305,0.6280] , D 2 = [0.6057,0.7859] , 

D 3 = [0.4979,0.7030] , D 4 = [0.6587,0.7894] . 

Step 4. Construct the likelihood matrix P using Eq. (11) and obtain the ranking vector 
co = (co 1 ,co 2 ,co 3 ,co 4 ) based on Eq. (12). 

0.5 0.0593 0.3231 0 

0.9407 0.5 0.7474 0.4091 

0.6769 0.2526 0.5 0.1317 ’ 

1 0.5909 0.8683 0.5 

rq =0.1569, a> 2 = 0.2998, co 3 - 0.2 134, and a> 4 = 0.3299. 

Step 5. According to the descending order of &>, (/ = 1,2, 3,4) , the ranking of all alternatives is 
a 4 >-a 2 ya 3 >- a t and a 4 is the best one. 

5.2 Comparison analysis and discussion 

In order to further verify the feasibility and effectiveness of the proposed approach based on INSs, a 
comparison analysis is conducted using three existing methods and the analysis is based on the same 
illustrative example. 

In Ref. [48] and Ref. [49], some aggregation operators were developed to aggregate interval 
neutrosophic formation. Two methods in Ref. [48] utilized the similarity measures based on the relationship 
with distances , where the similarity measure of Method 1 is on the basis of the Euclidean distance and the 
similarity measure of Method 2 is in view of the Hamming distance. The ranking results were got though 
calculating the similarity measures between each alternative and the PIS. The method in Ref. [49] relied on 
the Hamming distance between INSs and provided a compromise solution to rank alternatives, which 
considered the distance to both the absolute PIS and NIS. The comparison results can be found in Table 1. 



p= (p‘j L= 
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Table 1 



Comparison with different methods 



Methods 


Ranking vectors 






Ranking results 


Method 1 of Ref. [48] 


0.7370 


0.9323 


0.8344 


0.9034 


a 2 y a 4 y a 2 y a x 


Method 2 of Ref. [48] 


0.7667 


0.9542 


0.8625 


0.9600 


a 4 y a 2 y a 3 y a x 


Method of Ref. [49] 


0.941 


0.319 


0.666 


0.114 


a 4 y a 2 y a 2 y a x 


The proposed approach 


0.1569 


0.2998 


0.2134 


0.3299 


a 4 y a 2 y a 2 y a x 



The difference among the ranking results shown in Table 1 is the sequence of a 2 and a 4 . (1) The 

comparison methods in Ref. [48] were conducted using the similarity measures that only considered the 
relationship between each alternative and the PIS. There are some drawbacks. Firstly, it is not easy to choose 
the PIS in the assessment information given in the form of INSs. Secondly, the ideal solution is closely 
related to the number of alternatives as well as the evaluation values of alternatives. Thus it may vary as the 
original information changes. If only PIS is taken into account and NIS is ignored, the ranking of alternatives 
may be incorrectly reversed, and this may be amplified in the final results. There will be an identical ranking 
result with the proposed approach in case of replacing the PIS with the absolute one and allowing the usage 
of the absolute NIS at the same time, and utilizing the closeness coefficient (performance of each alternative) 
to determine the ranking of alternatives. The updated results are shown in Table 2. Therefore, the methods in 
Ref. [48] are not reliable enough, which can be demonstrated through the changes of the final results. (2) The 
result acquired by the fuzzy decision making method utilizing TOPSIS in Ref. [49] is consistent with the one 
obtained through Method 2 of Ref. [48] and the proposed approach, but the basis of these approaches is not 
similar. The method in Ref. [49] considered the distances to the absolute PIS and NIS based on the primary 
principle of TOPSIS. The performance of each alternative was aggregated in the form of a real number. 
However, unlike the method in Ref. [49], the proposed approach takes advantage of the merits of interval 
number and converts the performance of each alternative to an interval number, which can decrease 
information loss. Then these interval numbers are compared by establishing the relative likelihood-based 
comparison relations so as to rank all alternatives. What’s more, all methods in Ref. [48] and Ref. [49] were 
conducted by directly calculating similarity or distance measures between INSs, which are complex in 
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Table 2 



Ranking results by revising the methods in Ref. [48] 



Methods 


Ranking vectors 






Ranking results 


Method 1 of Ref. [48] 


0.5025 


0.6908 


0.5987 


0.6983 


a 4 y a 2 y a 3 y a x 


Method 2 of Ref. [48] 


0.5025 


0.6900 


0.5983 


0.6959 


a 4 y a 2 y a 2 y a x 



Furthermore, compared to the considered methods, the proposed approach makes good use of both 
interval numbers and TOPSIS method. It is well known that adopting interval numbers to represent 
performances in the TOPSIS method can give up various kinds of aggregation operators and score functions 
that may bring different results and effectively reduce the loss of decision information. Accordingly, the 
proposed approach is more flexible and reliable in handling MCDM problems than the compared methods in 
the interval neutrosophic environment. 

6. Conclusion 

INSs can flexibly express uncertain, imprecise, incomplete and inconsistent information that widely 
exist in scientific and engineering situations. So it is of great significance to study MCDM methods with 
INSs. In this paper, the basis definitions related to INSs are proposed and the useful operational laws and 
properties are discussed in detail. Then based on related research achievements in IFSs, a transformation 
operator and cross-entropy are defined. Thus, a MCDM method is established based on cross-entropy and 
TOPSIS, which computes the cross-entropy after transforming INNs into SNNs based on the transformation 
operator, and aggregates the performances of alternatives into interval numbers. Finally, the ranking result is 
obtained by comparing these interval numbers with a possibility degree method. 

The advantages of this study are that the approach can be both simple and convenient in computing and 
effective to decrease the loss of evaluation information. The feasibility and validity of the proposed approach 
have been verified through the illustrative example and comparison analysis. The comparison results 
demonstrate that the proposed approach can provide more reliable and precise outcomes than other methods. 
Therefore, this approach has much application potential in dealing with MCDM problems in the interval 
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neutrosophic environment, in which criterion values with respect to alternatives are evaluated by the form of 
INNs and the criterion weights are known information. 
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